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Предложена модель описания параметров состояния металлов, испытывающих мощное энергетическое воздействие. Показана 
работоспособность модели в широком диапазоне изменения плотностей и температур. Фундаментальность физических принци- 
пов положенных в основу модели, позволила получить новые результаты при рассмотрении классических задач гидродинамики. 


Развитие систем генерации сильноточных 
электронных пучков, мощных импульсных пучков 
ионов, лазерного излучения создало уникальную 
возможность получения новых видов воздействия 
концентрированных потоков энергии (КПЭ) на 
материалы. Процессы, сопровождающие взаимо- 
действие КПЭ с веществом, могут использоваться 
для решения большого числа научных и технологи- 
ческих задач, в том числе: получение сверхплотной 
плазмы инерциального термоядерного синтеза; 
модификация прочностных, триботехнических, 
антикоррозийных свойств металлов и сплавов; 
синтез новых композиционных материалов, вклю- 
чающих метастабильные фазы и соединения; полу- 
чение тонких пленок и покрытий посредством 
осаждения абляционной плазмы; разрушение и 
резка твердых материалов. Детально области при- 
менения КПЭ рассмотрены в монографии [1] и об- 
зоре [2], а также в библиографии к ним. 

Общим признаком всех областей использования 
КПЭ является возбуждение необходимого процесса 
в веществе, что реализуется определенным выбо- 
ром интенсивности и времени облучения, а также 
видом КПЭ и подбором оптимальных амплитудно- 
временных параметров импульса. Современные 
системы генерации КПЭ позволяют получать пото- 
ки в широком диапазоне интенсивности воздей- 
ствия Ж=10 5 ...10 13 Вт/см 2 . При длительности им- 
пульса КН*... КС 6 с материал может испытывать воз- 
действие с флюенсом КС... ІО 3 Дж/см 2 /на импульс. 
Такой широкий диапазон энергетического воздей- 
ствия КПЭ определяет возможность генерации в 
объеме мишени большого числа разнородных фи- 
зических явлений и, как следствие, многообразие 
результирующих эффектов, возникающих как при 
протекании отдельных процессов, так и при супер- 
позиции их совокупности. 


Используя в качестве критерия возникновение 
фазовых превращений в объеме мишени при воз- 
действии КПЭ, можно условно выделить три подди- 
апазона плотности мощности импульсного потока: 

1 . Для систем ионной имплантации и прямого нане- 
сения на поверхность наноразмерных пленок 
ионным пучком используются низкоэнергетичес- 
кие пучки (до десятков кэВ) и соответствующими 
плотностями тока, обеспечивающими плотность 
мощности до ІО 6 Вт/см 2 [1]. Выделившейся энер- 
гии, недостаточно для развития фазовых превра- 
щений, а модифицирующий эффект обусловлен 
формированием дефектной структуры в припове- 
рхностной области, не превышающей глубины 
проникновения ионов (десятки нм). 

2. Используемые в экспериментах по инерциально- 
му термоядерному синтезу КПЭ предельно высо- 
ких плотностей мощности (более ІО 11 Вт/см 2 ) при 
взаимодействии с веществом вызывают объем- 
ную сублимацию материала мишени. Рассмотре- 
ние динамики поведения мишени проводится в 
газовых и плазменных приближениях [2]. 

3. При рассмотрении КПЭ промежуточных плот- 
ностей мощности ІО 7 . ..ІО 10 Вт/см 2 в области 
энерговыделения и прилегающих к ней слоях 
мишени возможно формирование локальных 
зон различного фазового состава, в том числе с 
образованием слоев бинарных состояний "жид- 
кость-твердое тело" и "жидкость-пар". 

Наиболее многопараметрической в диапазоне 

промежуточных плотностей мощности является сис- 
тема "КПЭ-металл". На рис. 1 представлено характер- 
ное распределение зон различного фазового состава в 
приповерхностной области металлической мишени 
на момент окончания импульса ионного тока ( 10 7 с) 
и флюенсом около 30 Дж/см 2 . Несмотря на относи- 
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тельно малые геометрические размеры зон фазовых 
превращений, эти процессы оказывают существен- 
ное влияние, как на динамику формирования и раз- 
лета албяционного факела, так и на генерацию и эво- 
люцию импульса механических возмущений, расп- 
ространяющегося вглубь мишени. Последний может 
иметь сложную структуру (являться комбинацией уп- 
ругих пластических и ударных волн) и определять мо- 
дификацию свойств материала мишени далеко за 
пределами области теплового воздействия КПЭ (так 
называемый эффект дальнодействия [2]). 



Рис. 1. Распределение зон различного фазового состава. 
Области: 1) твердотельная; V) твердотельная, под- 
верженная пластическим деформациям; 2) газоплаз- 
менного облака, расширяющегося с облучаемой по- 
верхности; 3) переходных состояний; 4) расплава; 

5) бинарного состояния "жидкость -твердое тело"; 

6) бинарного состояния "жидкость -пар"; 7) возмож- 
ной модификации на тыльной поверхности при отра- 
жении импульса механической нагрузки 

Динамика процессов возбуждаемых КПЭ проме- 
жуточных плотностей мощности в настоящее время 
является наименее исследованной. Изучение всей со- 
вокупности явлений, сопровождающих такое воздей- 
ствие, представляет собой актуальную задачу не толь- 
ко технологического, но и фундаментального значе- 
ния. Применение методов механики сплошной среды 
позволяет рассмотреть поведение всей системы 
"КПЭ-металл" в целом. При этом, важное значение 
приобретает описание параметров состояния в широ- 
ком диапазоне фазовой диаграммы. Выбор коррект- 
ного широкодиапазонного уравнения состояния, яв- 
ляется необходимым для замыкания системы уравне- 
ний Навье-Стокса. Кроме того, для формирования 
законченной модели упруго-пластической среды не- 
обходимо рассмотрение полного напряженно-дефор- 
мированного состояния металла, его зарядового сос- 
тояния и детализации описания термодинамических 
параметров на границах раздела сред и фаз. 

Напряженно-деформированное состояние металла 

Описание динамики сплошной среды подразу- 
мевает формулирование законов сохранения с уче- 


том пространственно неоднородного поля механи- 
ческих напряжений, образованного интенсивным 
подводом энергии и деформацией вещества возни- 
кающей при движении среды. В каждой точке 
сплошной среды напряженное состояние характе- 
ризуется симметричным тензором напряжений, 
который можно разделить на две составляющие: 
шаровую часть и девиатор тензора напряжений 
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где о„ <7 У , сг, - нормальные, т ч =т ух , т ѵ =т а , т у =т, у - 
касательные напряжения. 

Шаровая часть тензора описывает гидростати- 
ческое давление р, которое изотропно. Девиатор 
тензора напряжений характеризует анизотропию 
напряженно деформированного состояния в об- 
ласти упругих деформаций, связанных с составля- 
ющими тензора напряжений законом Гука. 

При 1У< I О 6 Вт/см 2 изменение давления в облас- 
ти термализации пучка мало по сравнению с возни- 
кающими термическими напряжениями. В пренеб- 
режении гидростатическим давлением система 
уравнений Навье-Стокса преобразуется в волновое 
уравнение. При ІѴ>\(У Вт/см 2 прочностные свой- 
ства металла не являются определяющими, пренеб- 
регая девиатором тензора напряжений переходим к 
газодинамическому описанию динамики парамет- 
ров мишени. В промежуточном диапазоне плотнос- 
тей мощности необходимо учитывать обе составля- 
ющие тензора напряжений. При этом необходимо 
описать переход металла из упругого в пластическое 
состояние, например по условию текучести Мизеса 
(5 2 +5 , 2 2 +5 3 2 <(2/3)с 7 г 2 , здесь 5, - главные компоненты 
девиатора тензора напряжений; оу - динамический 
предел текучести). Это позволит получить правиль- 
ную интерпретацию свойств упругопластической 
среды вблизи состояния текучести, а именно полу- 
чить разрыв первой производной в критической 
точке на адиабатах Пуассона и Гюгонио. 

При описании термодинамических функций в 
широком диапазоне фазовой диаграммы, необхо- 
димо учитывать наличие сил межатомного взаимо- 
действия и их преодоление при расширении веще- 
ства, что в значительной степени определяет зако- 
номерности фазовых превращений. Соответствен- 
но давление разделим на две части: упругую состав- 
ляющую р„ связанную исключительно с межатом- 
ным взаимодействием, зависящим от расстояния 
между атомами; тепловую составляющую р т . В 
свою очередь тепловую составляющую будем раз- 
делять на слагаемые, соответствующие тепловому 
движению ионов р, и тепловому возбуждению 
электронов р е [3], 
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р(р,т) = рЛр) + рЛрН) = Ра (р) + Рі (рА)+р е (. р,т ), 
е(р, Г) = еДр) +е т {р,Т) = е 5 (р) +е, (р, Г) +е е (р, Г). 

Функциональные зависимости составляющих 
давления и соответствующих составляющих внут- 
ренней энергии от плотности и температуры опре- 
деляются уравнением состояния, выбор которого 
будет определяться следующими критериями: при- 
менимость для различных металлов и сплавов; ис- 
пользование только известных величин; минимум 
числа подгоночных параметров; термодинамичес- 
кая совместимость величин р иг. 


Уравнение состояния металла 

Ионные составляющие давления и удельной внут- 
ренней энергии определяются тепловым движением 
ионов, которое является классическим. В различ- 
ных агрегатных состояниях связь внутренней энер- 
гии с температурой определяется числом степеней 
свободы ионов, соответственно для: 

твердого тела е і = Зі^кТ 1 , ( де. 

г Рі = Р — 

идеального газа е . = (3 / 2)и ДГ ] др 



где щ=И а /М- концентрация ионов. Для выделения 
диапазона переходного состояния от твердого тела 
к идеальному газу использовался параметр 
Х = {Т т /Т) 1/ \ где Т т - температура Линдемана [4, 5]. 
В интервале 10^<^<1 обеспечивалась гладкая 
сшивка е,(р,7) с сохранением непрерывности пер- 
вой и второй производной на границах диапазона. 

Тепловое возбуждение электронов. Разделив теп- 
ловые компоненты давления и внутренней энергии 
на ионную и электронную составляющие, мы 
предположили, что металл представляет собой сво- 
бодный электронный газ, целиком заполняющий 
непрерывное пространство атомных ячеек, в цент- 
рах которых находятся ионы. При этом термодина- 
мические функции аддитивны, а соответственно, 
электронный газ рассматривается независимо от 
ионной компоненты. 

При 7=0 К электронный газ полностью вырож- 
ден. Электроны занимают самые низкие энергети- 
ческие состояния, а их кинетическая энергия не пре- 
вышает граничной энергии Ферми е Р Однако, эта 
кинетическая энергия включена в упругую энергию 
вещества, которая не имеет отношения к тепловой 
энергии [3]. При повышении температуры электро- 
ны частично переходят в более высокие энергетичес- 
кие состояния, что приводит к росту тепловой энер- 
гии электронного газа. В пределе высоких темпера- 
тур электроны перестают подчиняться статистике 
Ферми-Дирака, а электронный газ становиться иде- 
альным. Следуя авторам работы [4], введем параметр 


Пп е ,Т) 


Ъп е кТ 
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где л/р, 7)=2Др, 7)й,(р) - концентрация электро- 
нов, 2 т - средняя степень ионизации (количество 
электронов, приходящихся на одну атомную ячей- 
ку). Асимптотики параметра У определяют пре- 


дельные состояния электронного газа. Кроме того, 
важные термодинамические величины, такие как 
удельные свободная и внутренняя энергия могут 
быть выражены как функции только от п е и У 

Г = - Ѳ(У), е, =- ^Ч'(У), 

Р Р 


где ѳ(У) = |у 

р - химический потенциал электронного газа; 
Р х (а) - интегралы Ферми -Дирака. Параметр а не- 
явно зависит от У через уравнение 

(|г)"Ч,<«>-§. 

Асимптотика У— >со (идеальный газ) дает 
Т( У)~Ѳ( У)~ У, при У— >0 получим Ѵ Р(У)~-Ѳ(7)~1. В 
последнем случае для полностью вырожденного 
электронного газа Р— >е/р), тогда для свободной 
энергии электронного газа окончательно получаем 

4=^(р)- е ^[Ѳ(У) + 1]. 

Р 


а + 


2 ршісР 

3 Р ѵ1 (а) 


;ПЛ= 2 гЩ« = -Ь 

3 Р т (а) кТ 


Электронное давление и удельную тепловую 
энергию находим по свободной энергии из извест- 
ных термодинамических соотношений 


-И»! 


=р" 


(дР Л 

е_ 

а д р) т 


Упругие составляющие давления и удельной 
внутренней энергии зависят только от плотности и 
равны полному давлению и удельной внутренней 
энергии при абсолютном нуле температуры. В этом 
случае равновесное состояние твердого тела харак- 
теризуется взаимной компенсацией межатомных 
сил притяжения и отталкивания, а соответственно 
минимуму потенциальной упругой энергии, кото- 
рый можно принять за начало отсчета е=0 при 
этом плотность вещества р 0 , р= 0, Т=0. 

При расширении вещества преобладают силы 
притяжения. Силы взаимодействия быстро убывают 
при увеличении расстояния между атомами, поэто- 
му потенциальная энергия (упругая составляющая) 
возрастая, асимптотически стремится к постоянно- 
му значению равному энергии сублимации. При 
сжатии вещества преобладающую роль играют силы 
отталкивания, которые резко возрастают по мере 
сближения атомов, соответственно, е, быстро увели- 
чивается. Упругая составляющая давления связана с 
упругой составляющей энергии соотношением 


В литературе, например [4-7], рассмотрены 
различные методики расчета упругих компонент. 
При выборе таковой значимым является адекват- 
ное описание параметров состояния вблизи трой- 
ной точки. Именно в ее окрестности происходит 
разрыв межатомных связей, тепловая и упругая 
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составляющие давления (последняя отрицательна) 
примерно на порядок превышают результирующее. 
Для описания параметров состояния, использо- 
ванного нами, приемлемые результаты в области 
тройной точки дает интерполяционная зависи- 
мость, предложенная в [7] 
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где х=р 0 /р; а, р - параметры формулы; Д, - модуль 
объемной адиабатической сжимаемости. 

Ионизационный состав вещества в широком ди- 
апазоне фазовой диаграммы необходим для расче- 
та концентрации электронов не только при опреде- 
лении электронных составляющих термодинами- 
ческих функций, но и при описании процесса ре- 
лаксации энергии КПЗ во внутреннюю энергию 
среды. В металле даже при Т= 0 всегда присутствует 
некоторое количество свободных электронов. Во 
многих случаях для расчетов свойств разряженной 
многократно ионизованной плазмы хорошим 
приближением является модель "среднего иона" 
[3]. В [8] модель "среднего иона" обобщена на слу- 
чай больших плотностей и температур, в том числе 
для твердого тела. С учетом холодной ионизации 
вещества уравнение ионизации имеет вид 

/V {2 Я /Ѵ,Т) + 1{2 т )=Ь(2 /ѴУ(і + И Т 2 Ѵ°У\ 

где /и еГ (2„/Ѵ, Т) - химический потенциал идеального 
газа Ферми; К- объем атомной ячейки; І(І т ) - сгла- 
женная интерполяция дискретных потенциалов ио- 
низации; Ъ , Д /и, а - подгоночные параметры. 

При расчетах электронных компонент термоди- 
намических функций по ионизационному составу, 
полученному по данной методике, как и в [8], полу- 
чено существенное расхождение в значении давле- 
ния в окрестности тройной точки. Это обусловлено 
тем, что влияние ионизации давлением учитывает- 
ся дважды, один раз при определении электронных 
составляющих давления и внутренней энергии в 
концентрации электронов и второй раз косвенно 
при определении упругих составляющих. Поэтому 
для расчета параметров состояния использовалось 
эмпирическое соотношение, предложенное в [4] 


Т 

1 + 4^ 

г \ 

Р_ 

0,05 

Т 

<Ро ) 



где Г 0 = 4,32, эВ; р 0 =Д/ 10, г/см 3 ; 3,2- массовое чис- 
ло и заряд ядра атома соответственно. При опреде- 
лении параметров локальной области энерговыде- 
ления ионизационный состав определялся по 
обобщенной модели "среднего иона". 

Сравнение и обсуждение результатов расчета 

Результаты расчетов термодинамических пара- 
метров по предложенной модели сравнивались с 


расчетными и экспериментальными данными, 
представленными в [4-11] в различных диапазонах 
плотности и температуры для различных металлов. 
На рис. 2 сопоставлены результаты расчетов изо- 
терм меди, а также кривой плавления и линии рав- 
новесия жидкость пар с аналогичными данными из 
[9]. Кроме того, указаны экспериментальные пара- 
метры критической точки из обзора [11]. 


Р.кбар 



Рис. 2. Фазовая ( Р-р) диаграмма меди. Сплошные линии — 
предложенная модель; пунктир —данные из [9] 


При исследовании закономерностей процессов 
сопровождающих воздействие импульсных КПЗ 
на металлы, значительный интерес представляет 
корректное описании поведения вещества при пе- 
реходах через "критические" состояния. В частнос- 
ти, переход металла из упругого в пластическое 
состояние, формирование внутренних границ об- 
ластей фазовых переходов. На рис. 3 и 4 представ- 
лены результаты расчетов динамики параметров 
алюминиевой мишени, полученные в рамках лаг- 
ранжевой модели, реализующей схему Уилкинса 
[12], с детализацией начальных и граничных усло- 
вий [13] и предложенным в настоящей работе опи- 
санием параметров состояния. 

Рис. 3 иллюстрирует протекание ударно-волно- 
вых процессов в металле (алюминий). Рассмотрена 
модельная задача прохождения ударной волны (УВ) 
в упругопластической среде и ее отражение от сво- 
бодной границы. Закономерности этих процессов 
известны [3]. При прохождении по объему металла 
УВ с интенсивностью, превышающей предел теку- 
чести, наблюдается расщепление упругой и пласти- 
ческой волн, обусловленное разрывом первой про- 
изводной на адиабате Гюгонио в точке текучести и 
соответствующим изменением скорости звука. Вы- 
ход УВ на свободную границу сопровождается разг- 
рузкой вещества до нулевого давления и удвоением 
скорости вещества за фронтом отраженной волны. 

При моделировании прохождения произволь- 
ной УВ (рис. 3, а) по объему металла наблюдается 
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Рис. 3. Результаты моделирования выхода ударной волны на свободную границу в алюминии. Поля: а — полного напряжения ; 
р — плотности ; и — скорости 




Рис. 4. Формирование внутренней границы раздела фаз при воздействии на алюминиевую мишень мощного ионного пучка: 
а) динамика поля плотности в области термализации протонного пучка ; б) фазовые траектории лагранжевых ячеек, 
расположенных на границе парогазовой и конденсированной фаз (М — кривая плавления, Р — бинодаль) 
Примечание: символами О , Л , □ , * обозначены: а) положения лагранжевых ячеек и 6) их фазовое состояние в различ- 
ные моменты времени 


ожидаемая комбинация упругой и пластической 
волн (рис. 3 , б), при выходе которых на свободную 
границу впервые получен следующий результат. 


Сначала свободной границы достигает упругая 
волна, которая после отражения распространяясь 
вглубь мишени, встречается с пластической. Су- 
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перпозиция встречных пластической и отражен- 
ной упругой волн определяет формирование в пог- 
раничной области слоя вещества не подверженно- 
го пластической деформации (рис. 3, в). На рис. 3 
показано полное давление, равное сумме гидроста- 
тического давления и осевого компонента девиато- 
ра тензора напряжений, и, соответственно, равное 
полному напряжению, взятому с обратным знаком. 
Полное давление равно нулю за фронтом отражен- 
ной пластической волны, а значение скорости со- 
ответствует правилу удвоения. Необходимо отме- 
тить, что в распределении плотности вещества за 
фронтом отраженных волн, наблюдается одноз- 
начное разделение областей упругой и пластичес- 
кой деформации. В области пластического течения 
формируется остаточная деформация [6]. 

Для моделирования процесса испарения металла 
при мощном энергетическом воздействии особое 
внимание уделяется вопросу о границах раздела фаз. 
Например, в [5], в качестве такой границы выбрана 
бинодаль и предложены специальные граничные ус- 
ловия для определения значений удельного объема 
(плотности) и равновесного давления. На рис. 4, а 
представлена динамика поля плотности в области 
локального энерговыделения, обусловленного тер- 
мализацией мощного протонного пучка (амплитуд- 
но-временные параметры пучка соответствуют [14], 
при Ѵ тш = 660 кэВ,у шах =1,5 кА/см 2 , г=120 нс) в алю- 
миниевой мишени. К моменту окончания импульса 
ионного тока на глубине порядка 6 мкм формирует- 
ся локальная граница между парогазовой и конден- 
сированной фазами, характеризуемая разрывом в 
функции р(і). На рис. 4, б, на фазовой диаграмме 
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представлены траектории соседних лагранжевых 
ячеек, расположенных в окрестности разрыва. Вид- 
но, что область бинодали является внутренней гра- 
ницей, образовавшейся при разрыве межатомных 
связей на фазовом переходе жидкость-пар. Конден- 
сированное вещество в окрестности новой границы 
начинает двигаться вглубь мишени и испытывает 
разгрузку. Оторвавшееся парогазовое облако испы- 
тывает адиабатическое расширение, подчиняюще- 
еся законам газовой динамики. 

Заключение 

Предложенная модель адекватно и термодина- 
мически полно описывает параметры состояния 
металла в широкой области фазовой диаграммы. 
Приведенное сравнение результатов показывает 
приемлемую точность расчетов для задач модели- 
рования динамики сплошной среды. Корректное 
описание напряженно-деформированного состоя- 
ния позволило получить новые результаты. В част- 
ности, при рассмотрении задачи о выходе последо- 
вательности упругой и пластической волн на сво- 
бодную границу обнаружено формирование в пог- 
раничном слое области вещества, не подверженной 
пластической деформации. Использование модели 
в расчетах динамики системы "КПЭ-металл" дает 
возможность описания процессов фазовых превра- 
щений, в том числе однозначного определения гра- 
ниц областей различного фазового состава. Пос- 
леднее не требует постановки специальных усло- 
вий на внутренних границах раздела фаз и является 
свойством разработанной модели. 


9. Альтшулер Л. В., Бушман А.В., Фортов В.Е. и др. Изэнтропы 
разгрузки и уравнение состояния металлов при высоких плот- 
ностях энергии //ЖЭТФ. — 1980. — Т. 78. — № 2. — С. 741—760. 

10. Замышляев Б.В., Ступицкий Е.Л., Гузь А.Г., Жуков В.Н. Состав 
и термодинамические функции плазмы: Справочник. — М.: 
Энергоатомиздат, 1984. — 144 с. 

1 1 . Басин А.С. Основные параметры критической точки металлов 
с плотноупакованной кристаллической структурой // Химия и 
компьютерное моделирование. Бутлеровские сообщения. — 
2002. — Приложение к выпуску № 10. — С. 103—108. 

12. Уилкинс М.Л. Расчет упруго-пластических течений / В сб.: 
Вычислительные методы в гидродинамике. — М.: Мир, 1964. — 
С. 212-263. 

13. Бойко В.И., Данейкин Ю.В., Юшицин К.В. Численное описа- 
ние процессов генерации волнового возмущения, обусловлен- 
ного воздействием на металлическую мишень импульсного 
энергетического потока // Вопросы атомной науки и техники, 
серия "Физика радиационных повреждений и радиационное 
материаловедение". —2003. — № 3. — С. 133—137. 

14. Бойко В. И., Быстрицкий В.М., Волков С.Н. и др. Генерация и 
фокусировка мощного ионного пучка в магнитоизолирован- 
ном диоде // Физика плазмы. — 1989. — Т. 15. — Вып. 11. — 
С. 1337-1341. 


21 




